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ABSTRACT 

Studies based on high-resolution spectroscopic data of rapidly oscillating Ap stars 
show a surprising diversity of pulsation behavior in the atmospheric layers, pointing, 
in particular, to the co-existence of running and standing waves. The correct inter- 
pretation of these data requires a careful modelling of pulsations in these magnetic 
stars. In light of this, in this work we present a theoretical analysis of pulsations in 
roAp stars, taking into account the direct influence of the magnetic field. We derive 
approximate analytical solutions for the displacement components parallel and per- 
pendicular to the direction of the magnetic field, that are appropriate to the outermost 
layer. From these, we determine the expression for the theoretical radial velocity for 
an observer at a general position, and compute the corresponding pulsation amplitude 
and phase as function of height in the atmosphere. We show that the integral for the 
radial velocity has contributions from three different types of wave solutions, namely, 
running waves, evanescent waves, and standing waves of nearly constant amplitude. 
We then consider a number of case studies to illustrate the origin of the different pul- 
sational behaviour that is found in the observations. Concerning pulsation amplitude, 
we find that it generally increases with atmospheric height. Pulsation phase, however, 
shows a diversity of behaviours, including phases that are constant, increasing, or de- 
creasing with atmospheric height. Finally, we show that there are situations in which 
the pulsation amplitude goes through a zero, accompanied by a phase jumps of tt, and 
argue that such behaviour does not correspond to a pulsation node in the outermost 
layers of the star, but rather to a visual effect, resulting from the observers inability 
to resolve the stellar surface. 
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1 INTRODUCTION 

The rapidly oscillating Ap (roAp) stars are found among 
the coolest subgroup of chemically peculiar Ap stars, which 
is located in the main sequence part of the classical insta- 
bility strip. They pulsate with p eriods typica l ly wit hin the 
range from 5 to 21 minutes, (e.g. iKurtz et all (|2005l )). have 
oscillations with amplitudes between 0.5 and 5 kms -1 in 
velocity, and, as in other Ap stars, have strong, large scale 
magnetic fields, with typical intensities of a few kG, although 
in some st ars the magnetic field strength can be higher than 



20 kG.fe.g. lHubrig et al 
two solar masses (|Kurt j 



(2009)). Moreover, they have about 



1990) , and temperatur es that range 



from about 6400 to 8100 K (|Kochukhovl (|200gh ). 

The present number of known roAp stars is more 
than 40. Due to their characteristics, with roAp stars 



we have the unique opportunity to observe the interac- 
tion of acoustic modes with strong large-scale, magnetic 
fields. The roAp stars ha ve been obser ved photometrically 
since their discovery by iKurt j (JI982) and until recently 
the knowledge about their acoustic oscillations was essen- 
tially based on high-speed, ground-based, photometric ob- 
servations. However, in the past few years numerous excit- 
ing observational results have been published, both as a 
result of the acqu i sition of space-based ph o tometry (e.g. 



Gruberbauer et all 120081 ; iBruntt et all 2009; B alona et all 



2010T) and as a resu lt of the analysis of h igh-resolution spec- 



trosc o pic data (e.g. Kurtz et a l. 2006a b; iRvabchikova et all 



120071 ; Isachkov et alj 120071 )7 Such high-resolution spectro- 
scopic data hold unique information about the structure 
and dynamics of the peculiar atmospheres of roAp stars, 
and have revealed a surprising diversity in the pulsation be- 
havior of different lines in the roAp spectra. 



* E-mail: jsousa@astro.up.pt 



The general picture that emerges from the analysis of 
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time-series of high-resolution spectroscopic data of roAp 
stars is that pulsational variability is seen predominantly 
in lines of rare-earth ions, especially those of Pr and Nd, 
which are strong and numerous in the roAp spectra and 
are formed in the higher layers of the atmosphere. On the 
other hand, lines of light and iron-peak elements, enhanced 
in the lower atmospheric layers, often do not show pulsation 
variability withi n the observational detect ion limit for most 
roAp stars (e.g. iRva bchikov a et all 12002) ') . Thus, generally 
pulsations are found to be weak or non-detectable in the 
lower atmosphere, while often reach amplitudes of several 
kms -1 higher in the atmosphere. 

Another interesting feature that can be inferred from 
spectroscopic observations of roAp stars is the presence of 
significant shifts in pulsation phase when comparing ra- 
dial velocities derived from lines of different rare-earth ele- 
ments, or even from different line s of the same element (e.g. 
iKochukhov fc Rvab chikova 200j]). Likewise, phase shifts are 
seen when performing d epth- in- line analys is of the Ha and 
Nd III 6145 A lines, (e.g. lKurtz et al.ll2006al ). Along with the 
phase shifts, depth variations of pulsation amplitude are also 
inferred, which in the generality of the cases point to an in- 
crease of amplitude with atmospheric height. Moreover, in 
some stars the radial velocity amplitudes measured in lines 
of the rare-earth ion Pr III and in the core of the Ha line, 
are seen to vary on a time scale of a few pulsation cycles 
l|Kurtz et al.ll2006bl ). 

The phase behaviour derived from spectroscopic stud- 
ies of roAp stars is usually interpreted as resulting from 
the presence of running waves, or standing waves, in the 
atmosphere, depending, respectively, on whether, or not, 
the phase sh ows a variation with dept h. In particular, it 
was fo und by ISachkov et al.l (|2007l ) and IRvabchikova et all 
l|2007l ). that for stars with pulsation frequencies below the 
acoustic cut-off frequency, the pulsations seem to have 
standing wave character in the deeper layers and then, in 
the outer layers, behave like running waves. On the other 
hand, for stars which have pulsation frequency close to, or 
higher than, the cut-off frequency, the authors found that 
the pulsations behave like running waves from the deepest 
layers. 

One important and interesting feature that can be seen 
in spectroscopic data of roAp stars is the existence of pulsa- 
tion zeros in the amplitude diagram accompanied by phase 
jumps in the phase diagram. This was interpreted as the 
presence of a pulsation node at some depth in the atmo- 
sphere. T he first alleged node found in a roAp star was de- 
tected by iBaldrv et al.l (| 19981 ) in the atmosphe re of a Cir. 
More recently, in the case of the star lOAql, lElkin et al.l 
(2008) found that the pulsation phase derived from the lines 
of Tb III and Dy III differ by 180 ° . According to the authors 
this could be an indication of a radial node between the line- 
forming layers of these elements. Moreover, they found that 
the line bisectors for strong Ndlll line profiles show signifi- 
cant changes of phase, or even phase jumps, as in the case of 
the line Ndlll at 5102.41 A. For this line, the authors found a 
phase shift of 180° and at the same depth, an amplitude close 
to zero , which suggests the exis tence of a radial node. In ad- 
dition, iMkrtichian et all |2003l ) studied the radial velocity 
variations in the roAp star 33 Librae. The authors found 
that the Ndlll lines pulsate nearly in anti-phase with those 
of Ndll, a feature that the authors attribute to the presence 



of a pulsation node in the atmosphere of this star. The pul- 
sational behaviour of t he Nd ii and Nd iii lines in 33 Lib has 
also been discussed bv lKurtz et"ail <|2005l) . They found that 
the line depth versus pulsation amplitude supported the hy- 
pothesis of a node between t he lin e-forming layers of these 
two ions. IRvabchikova et al.l (|2007l ) found similar results for 
this star, confirming a phase jump between the Ndlll and 
Ndll lines, and the amplitude decreasing towards zero at the 
same atmospheric height. 

The results of the our theoretical analysis will be pre- 
sented in two separate papers, as detailed in sec. [2] In sec. [3] 
we discuss the physics underlying our model, along with the 
assumptions made. In sec.|?]we derive the expression for the 
theoretical radial velocity and describe a Toy Model that 
will be useful for the interpretation of the results. The re- 
sults of the analysis of six case studies are presented in sec.[S] 
followed by a general discussion in sec. [6] 



2 PARAMETER SPACE 

Given an underlying equilibrium stellar model, the general 
problem that we set out to study depends on a number of 
input parameters. Some of these are intrinsic to the star, 
namely: the magnetic field intensity and topology, charac- 
terized by the vector B; the unperturbed (i.e., in the absence 
of a magnetic field) oscillation mode, characterized by its ra- 
dial order, n, angular degree, I, and the azimuthal order m; 
the location of the chemical elements whose spectral lines 
are used in the derivation of the radial velocity, character- 
ized by an atmospheric depth (or region in depth) , as well as 
by longitudinal and latitudinal limits. Moreover, some pa- 
rameters depend on the position of the observer, such as: 
the inclination angle between the magnetic axis and the di- 
rection of the observer 8 (at a given phase of rotation); the 
inclination angle between the latter and the rotation axis, i. 

In our analysis we will keep some of the above fixed. 
In particular, we will only consider magnetic fields of a 
dipolar topology. We will not consider the effect of rota- 
tion on the dynamics and, consequently, will assume that 
the magnetic and pulsation axis are aligned. This is gen- 
erally a good approximation for roAp stars, except possi- 
bly in c ases when the manetic field is rather low (below 
«1 kG) l|Bigot fc Dziembowskill2002l ). Moreover, the unper- 
turbed pulsation modes to be considered are only those ax- 
isymmetric about the magnetic axis, thus, characterized by 
m — in a spherical harmonic decomposition about that 
axis of symmetry. Finally, we will consider only a fixed phase 
of rotation, so that our problem will have no dependence on 
the inclination angle i. We note that the relaxation of most of 
the above conditions is relatively straightforward, although 
in some cases will lead to significant additional work. 

The results of our study will be presented in two papers. 
The present paper will deal with the underlying mathemat- 
ical analysis and the in-depth investigation of a number of 
test cases - thecase studies discussed in section [5] For that, 
we will restrict the parameter space further, by considering 
only modes of degree I = 1 and the case in which the pulsa- 
tion and magnetic axis is along the line-of-sight (i.e., S — 0). 
The latter is justified by the fact that in this case the spheri- 
cal coordinate system associated with the observer coincides 
with the one used to solve the pulsation problem in the star. 
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Figure 1. Schematic representation of the star, vertical cut. 



As will become clear in sec. 14.21 that is a necessary condition 
for the in-depth study that we aim at in the present paper. 

In a second paper we will present the conclusions of 
a systematic search of the parameter space, including also 
variation of mode degree and of the position of the observer. 
Additionally, in that paper we will compare our results with 
those derived by simulations of line profiles, as well as with 
general trends of pulsation phase and amplitude found in 
the observations. 



3 PULSATIONS IN THE PRESENCE OF A 
MAGNETIC FIELD 

The study of pulsations in the presence of a magnetic 
field has been undertaken by a number of authors in the 
past, adopti ng different approximations , or mathematical 
approaches IIDz icmbowski & Goodc 1996; B igot et al.l | 2000 |; 



Cunha fe Goughl l200d ; ISaio fc Gautschvl 12004 ISaiol 120051 ; 
Cunhall2006l ) 

The magnetic field has a direct influence on the oscilla- 
tions only in the region where the Lorentz forces are com- 
parable to or larger than the gas pressure gradient, that is 
to say, in the surface layers of the stars. Therefore, we fol- 
low the approach of previous studies, and divide the star 
into two regions, a region where the magnetic pressure is 
greater than or comparable to the gas pressure, the surface 
boundary layer, and a region where the gas pressure is much 
larger than the magnetic pressure, the interior, defined to be 
the region between the center of the star and the boundary 
layer. In the interior the dynamics is not directly affected by 
the magnetic field and the oscillations there will be modified 
only as a consequence of the change in the conditions of the 
surface layers. 

The surface boundary layer itself can be divided into 
two regions, one where the magnetic pressure is comparable 
to the gas pressure, that we shall call the coupling region, 
and one where the magnetic pressure is much larger than the 
gas pressure, which is the magnetically dominated region. 

Concerning the oscillations, the studies referred to 
above have shown that in the coupling region the magnetic 



and acoustic wave components are coupled into a magnetoa- 
coustic wave. Moreover, the authors have shown that in the 
interior of the star, where the gas pressure is much larger 
than the magnetic pressure, the magnetoacustic wave de- 
couples into a fast (acoustic in nature) and a slow (magnetic 
in nature) component, the latter being an inwardly propa- 
gating wave. Likewise, higher in the atmosphere, where the 
magnetic pressure is much larger than the gas pressure, the 
magnetocaoustic wave decouples into acoustic and magnetic 
components, the lat ter being a standing w ave with almost 
constant amplitude |Sousa fc Cunhdl2008h . 

In Fig. [T] we show a schematic representation of a verti- 
cal cut in the star. On the left hand side we show the interior 
and the surface boundary layer and on the right hand side 
we show a zoom of the surface boundary layer in which the 
coupling and magnetically dominated regions are depicted. 
Moreover, the figure shows also the outer atmospheric re- 
gion, which we define as a subsection of the magnetically 
dominated region in which the plasma is assumed to be ap- 
proximately isothermal. In this study we will concentrate 
in the outer atmospheric region, and will take as input the 
velocity field at the base of that region which, in turn, is 
obtained from the study of t he star as a whole, following 
the approach of lCunhal ((2006). 



3.1 Surface boundary layer 

In our calculations we follow ICunha fc Goughl (2000) and 
assume a dipolar magnetic field, B, of polar magnitude B p 
which is force free, and hence does not influence the equilib- 
rium state. Therefore, in this layer the equilibrium state is 
governed by the following system of equations, 



dpo 

dt 
dp 

dt 

Vp = Po9o, 



0. 



0, 



(1) 

(2) 
(3) 



where p is the pressure, p is the density, g is the gravitational 
field, t is the time, and the subscript denotes the equilib- 
rium quantities. The equilibrium gravitational field, g , can 
be written as the gradient of the gravitational potential, <f>, 



go 



-v$, 



where $ satisfies Poisson's equation, 
V 2 $ = A-nGpo, 



(4) 



(5) 



where G is the gravitational constant. 

In the limit of perfect conductivity, adiabatic motions 
associated with a velocity, v, are governed by the following 
system of magneto-hydrodynamic equations, 



dB 

dt 
dp 



= V x (v x B), 
- pV.v = 0, 



dt 

dp _ 7p dp 



dt 



p dt ' 



dv _ . „ 

p— = -Vp + j xB + pg, 



(6) 
(7) 
(8) 
(9) 
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where d stands for total derivatives, d stands for partial 
derivatives, 7 is the first adiabatic exponent and j is the 
current density. 

Assuming small adiabatic perturbations of the equilib- 
rium state, such that, 

? = Qo + qi, 

for any scalar or vector component, q, and neglecting the 
terms of higher order in the perturbations, the two previous 
systems of equations can be rearranged to give a new system 
of equations which govern linear adiabatic pulsations in the 
Cowling approximation, namely, 



<9 2 £ 

Po-gjj = -Vpi + pigo + — (V x Bi) x Bo, 



1 



Pi 



-£.Vp - 7PoV.£, 



pi = -f.Vpo - poV.£, 
B x = V x (fx Bo), 



(10) 

(11) 
(12) 
(13) 



where po is the magnetic permeability and £ is the vector 
displacement, defined by the relation vi=d£/dt, and where 
we used the fact that the magnetic field in the equilibrium 
state is irrotational. 

As expected, the magnetic term in equation (|10[) is per- 
pendicular to the direction of Bo and, thus, within the ap- 
proximations considered, there is no magnetic component 
of the restoring force along the direction of the unper- 
turbed magnetic field. To proceed we follow the assump- 
tions adopted in previous studies and consider, at each 
latitude, a local plane-parallel layer with locally constant 
Bo = [B x ,0,B z ] and g , and neglect dBo/dx and dBo/dz, 
but allow B to change with latitude. In the above, we con- 
sidered a local coordinate system (x,y,z), with x and y point- 
ing along the latitudinal and longitudinal directions, respec- 
tively. Moreover, 2 is taken to be zero at some fiducial radius 
(that shall be identified later) located in the outer atmo- 
spheric region, and increasing outwardly. 



3.1.1 Components of the displacement 

Let us first begin with the projections of equation (|10|) in the 
directions parallel and perpendicular to the magnetic field. 
For that we will define £|| and £± such that £ = £^e^+£±e±, 
with en = Bo/|Bo|, and e± — (Bo x e y )/\Bo\, the unit 
vectors in the direction parallel and perpendicular to Bo, 
respectively, and e x and e y the unit vectors in the direction 
of the x and y axes, of the adopted Cartesian coordinate 
system, respectively. We note that due to the axisymmetry 
of both the magnetic field and the oscillation modes to be 
considered, the problem reduces to fourth-order (two second 
order coupled differential equations), and the only projec- 
tions that need to be considered in the analysis lie in the 
(e x , e z ) plane. 

Since the magnetic term in equation (|1UI) is zero along 
the direction of the magnetic field, when combined with 
equations (|11|) and (|12[) . the projection of equation (|10l) in 
this direction becomes, 



Po 



dt 2 



= 7V||(poV.£)+ J o <?(e !c .V)6 



(14) 



where V|| is the directional derivative in the direction par- 
allel to the magnetic field which is denned as V| = en.V, 
and go = ge z . 

Next, we assume that all horizontal derivatives of the 
displacement, and its components, are much smaller than 
the corresponding vertical derivatives. For low degree, high 
radial order modes in the absence of a magnetic field this is 
clearly a good approximation. In the presence of a magnetic 
field there are additional dependencies of the cigcnfunctions 
on the latitude, but these occur on a scale of 1/R, where R 
is the radius of the star, and thus introduce local horizontal 
derivatives that are much smaller than the vertical deriva- 
tives in the atmosphere. Thus, this approximation remains 
good in general, in the presence of a magnetic field. Note, 
however, that this approximation may fail, at particular lat- 
itudes, for particular frequencies, when the phase and/or 
amplitude varia tions of the displaceme nt with latitude are 
very sharp (e.g. lCunha fc Goughl koOO). However, since the 
observer sees only an integral of the eigenfunctions over the 
stellar disk, that failure is likely not to have a strong impact 
on the results, unless the eigenfunction is trapped within a 
given latitudinal region, in which case the approach consid- 
ered here will not be adequate. 

So, neglecting the horizontal derivatives of the displace- 
ment when compared with the vertical derivatives, we get, 
after some algebra, the following expression for the projec- 
tion of equation (|10[) in the direction along the magnetic 
field, 



Po|B | 2 d 2 C|| 
7Po-B| dt 2 



<9 2 £|| B x <9 2 C± 



dz 2 



B z dz 2 



po dz \ dz B z dz 



(15) 



Let us now consider the projection of equation JTU} in 
the direction perpendicular to the magnetic field. Taking the 
scalar product with e± we get, 



Po 



dt 2 



-Vpi + pigo] 



— |B | V||Bi± 

p>o 



\Bo\ VxBn 



(16) 



where Vi is the directional derivative in the direction per- 
pendicular to the magnetic field, defined as Vj_ = ex-V, 
and Bix and Bi| are, respectively, the components of Bi 
along and perpendicular to the magnetic field direction. 

In equation (|16[1 . the first term on the right hand side 
corresponds to the restoring force that would act in the 
absence of the magnetic field while the second term cor- 
responds to the magnetic response. After developing both 
terms and combining with equations (|11[) . (|12[) and (|13[) . 
equation (|16l) becomes, 



po- 



9 2 a |B | 2 2 



dt 2 im> 



V 2 £± +7Vx(poV.£)+ 



„ r B x _ -z 

Pod 1 — VxCn - — V||£|| 



(17) 



Bo I I Bo I 

Furthermore, similarly to what we have done for the 
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projection in the direction along the magnetic field, we ne- 
glect the horizontal derivatives of the displacement com- 
ponents when compared with the vertical derivatives, and 
equation (|17|) becomes, 



popo d 2 £± 
I Bo 1 2 dt 2 



1 + 



B~ 



d 2 L 



I Bo I 



B X B Z (3 dp 
I B 1 2 Po dz 



I Bo I 



dz 2 



dz B 



dz 



(18) 



where j3 = 7p Mo/|B | 2 . 

We can combine equations (|15[) and (jTHJ) and rewrite 
the system as follows, 



&% B x d 2 l 



dp d£\\ 



+ 



B z dz 2 podz dz 
d Po B x di± |B () | 2 poS 2 C|| 



podz B z dz B 2 -ypo dt 2 



a Po d 2 U 

dz 2 jpo dt 2 



7Po B z dt 2 



(19) 



(20) 



So, together, equations (|19p and (|20|) govern the dis- 
placement components parallel and perpendicular to the di- 
rection of the magnetic field, derived in the plane-parallel 
local analysis and under the approximations presented in 
this section. 



3.2 Magnetically dominated region 

There is a close relation between the location of the cou- 
pling region in a roAp star and the strength of the magnetic 
field in that star. The photosphere of a typical roAp star is 
located in different parts of the surface boundary layer, de- 
pending on the magnetic field intensity. For magnetic fields 
of the order of 1 kG, it is located in the coupling region, 
while for magnetic fields of the order of 3 kG it is located 
in the magnetically dominated region !|Cunhall2007h . There- 
fore, the oscillations in the atmospheres of these stars might 
look significantly different from star to star, and in different 
regions within the same star. In our analysis we are partic- 
ularly interested in the form of the velocity in the outer at- 
mospheric region. This means that in the region of interest, 
the magnetic pressure is much larger than the gas pressure 
and, consequently, the magnetoacoustic wave is decoupled 
into its acoustic and magnetic components. 

Thus, we next consider regions where is much smaller 
than one, and assume the displacement can be expressed as 
the sum of a slow (essentially acoustic) component, £ s , and 
a fast (essentially magnetic) component, £/ , i.e., 



(21) 



such that (f si ) _1 df si /dz >> (£f i )~ 1 d£ fi /dz, where £ si and 
£/i stand, respectively, for any of the components of the vec- 
tors £ s and £/. If the two components are indeed decoupled, 
equations (|19p - (|20p must be satisfied separately by £ s and 
£f. Moreover, we will consider solutions of the form, 



i(x,z,i) =£(x,z)e 



where to is the oscillation frequency, so that = — oj 2 £. 
Finally, we note that, as discussed in the introduction, spec- 
tral lines of the rare earth elements form very high in the 
atmosphere where the it becomes close to isothermal. Since 
our ultimate aim is to compare the velocity field derived 
from this analysis with that deduced observationally from 
the lines of the rare earth elements, we will make an ad- 
ditional approximation and consider only those outermost 
layers where we may expect the isothermal approximation 
to be appropriate. Also, we take the first adiabatic exponent, 
7, to be constant in that region, and equal to 5/3. 

In these layers we can defin e the pressure, po, and the 
density, po, as (e.g. iGoughll 19931) 



Po = Poe 



po = poe 



z/H 



-z/H 



(22) 
(23) 



where po and po are the pressure and density, respectively, at 
the base of the isothermal region, and H is the scale height 
defined by H^ 1 = — 1/podpo/dz = —1/ podpo/dz. With these 
expressions we simultaneously define z = as being the 
position of the base of the isothermal region. 



3.2.1 Slow Component 

In the limit when the fast and slow components are decou- 
pled, the approximate form of the space-dependent part of 
the slow component, £ a = £i| 4 ei| + £,± s e± can be obtained 
from the system of equations, 

d 2 £|| s dpo d£|| s bJ 2 p |B | 



dz 2 

d 2 ^ 
dz 2 
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Podz dz 
w 2 po Bi 



7Po Bj 



IPO 



£ik = o. 



B 2 



-f||s=0. 



(24) 
(25) 



These are derived from the system of equations (|19[) - (|20|) . 
after neglecting small terms (as discussed in Appendix A). 
Moreover, the latter implies that << £|| S) as expected 
for a wave that is essentially acoustic in nature, in a mag- 
netically dominated region. 

In the present case, of an isothermal atmosphere, equa- 
tion (|24l) admits solutions of the type, 



t Ua = Po - 1/2 (Ae^ +Be~ 



where, 



Ml 



1 dp \ 2 1 d 2 p 

2 podz J 2po dz 2 



u) 2 p |B | 2 



IPoB 2 



(26) 



(27) 



and A and B axe complex, depth independent, amplitudes. 

From equation (127f> we can derive the acoustic critical 
frequency, u ac , 
We find, 



defined as the frequency at which kf, — 



— cos(a z 



(28) 



where a z is the angle between the direction of the local mag- 
netic field and the local vertical coordinate. If we consider 
the case of the magnetic pole (a z = 0), we find that the 
latter equation is similar to the cutoff frequency for acoustic 
waves in the absence of the magnetic field. Therefore, here- 
after we shall use acoustic cutoff frequency for the acoustic 
critical frequency at the magnetic pole, and retain the name 
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acoustic critical frequency for the critical frequency at a gen- 
eral latitude. 

Going back to the solution expressed by equation (|26[) , 
we find that when the oscillation frequency is larger than 
the acoustic critical frequency, A must be set to zero, as 
otherwise energy would be sent in from outside the star. In 
that case, we find, from equations 1)251) and (|26l) . 



f||.= 



e H z 

1/2 
Po 



6 



_ ~ uj 2 p B. 



Be 



7p 3/2 B z (<fc|, + l/(2£0) a 



(29) 



If we consider the case of an oscillation frequency smaller 
than the acoustic critical frequency, we find that B must be 
zero, so that the energy does not increase with z. In that 
case, equations (|25jl and (|26|) imply that, 



T/2 e 
Po 



6 



pU 2 po B. 



Ae l 



7p 3/2 Bz + l/(2ff)) s 



(30) 



In both cases, we have set the integration constants that 
result from the integration of equation (|25|) to zero, as 
non-zero constants would correspond to solutions that do 
not obey the condition set for the slow wave, namely that 

(Zj^dU/dz » (tfi^dtfi/dZ. 

From equation (I28p we can see that oj ac depends on co- 
latitude, because it depends on the inclination of the mag- 
netic field. So, even when the oscillation frequency is lower 
than the acoustic cutoff frequency, there will be a colatitude 
at which, due to the inclination of the magnetic field, the 
acoustic wave will start to propagate. For this colatitude, 
above which the acoustic component starts to propagate as 
a running wave, we define a critical angle, 9 ac . 



3.2.2 Fast Component 

Similarly, in the limit when the fast and slow components 
are decoupled, the approximate form of the space-dependent 
part of the fast component, £/ = £||/e| + £±/ej_ can be 
obtained from the system of equations, 



d 2 £,±i , dp d£± f \B \ 2 uj 2 p ^ 

-I ; ; 1- - ?n/ - u - 



dz 2 podz dz 

-^- + *x€±/-o, 



B x B z jp 



where, 



2 

u po 
ypo 



(31) 
(32) 

(33) 



and va is the Alfven velocity. 

Equations (|31|l and (|32[) are derived from the system 
of equations (|19[) - (|20[ ). by neglecting small terms (see Ap- 
pendix A for details). An additional conclusion of that anal- 
ysis is that << £j_f, as expected for a wave that is 
essentially magnetic in nature, in a magnetically dominated 
region. Moreover, we note that k\ is always positive, which 
means that unlike the case of the slow component, in this 
case there is no magnetic critical frequency. 

In the present case, of an isothermal atmosphere, equa- 
tion (|32|l admits solutions of the type, 



C-L/ = Ci Jo I 



/ 2Huj^Jm) 1/2 



I -Bo 



Po 



+ C 2 Y Q 



( 2Hu) y fpo~ . 
{-lBoT P ° 



-(34) 



where Ci and C 2 are depth independent amplitudes. Since 
the function Yo diverges when its argument tends to zero, 
Yo diverges as po — > and therefore Ci must be zero to 
obtain a physically meaningful result. Moreover, when j3 is 
sufficiently small, equation (|31[) can be further approximated 
by, 



1 dp rig 

po dz dz 



— 1 ppU) po 

ft B x Bz 



(35) 



Thus, we find that the components of the fast displacement 
are, within the approximations considered, 



Ci Jo 



Ci- 



/2Hu 



\Bo\ 

PB X B. 



MO 1/2 

Po 



H \B \uip 



1/2 



( 2Hu)y/p7; 1/2 

Jl p I Po 
po V l-fol 



(36) 



Since the density po tends to zero, as one considers 
higher layers in the atmosphere, Jo tends to one there. Con- 
sequently, the solution for the fast component of the dis- 
placement perpendicular to the direction of the magnetic 
field, will tend to a constant. Moreover, both /3 and Ji tend 
to zero, as one considers higher layers in the atmosphere. 
Thus, the solution for the fast component of the displace- 
ment parallel to the direction of the magnetic field will tend 
to zero, as it should, since the pure magnetic wave will have 
a displacement perpendicular to Bo. 

3.3 Dimensionless Equations 

Prior to computing the expected radial velocity, we convert 
the quantities derived to this point that will enter such com- 
putation to corresponding dimensionless quantities. To that 
end we define: 

n = z/R, (t = oj/ljo, p = po/Po, p = po/po, 

£ N =^/R, s_l=U/R, C s =p /(p uoR 2 ),- 

h = B./ipop^lR 2 ) 1 ' 2 , b = \B \/( popoujlR 2 )^ 2 , 
where R is the radius of the star and ujo = V G 'M / R? , M is 

the mass of the star and G is the gravitational constant. We 
note that p and p defined above are not to be confused with 
the dimensional quantities named by the same symbol and 
used in equations (JB])-®. With these definitions the dimen- 
sionless scale height becomes W~ 1 = R/H = ~p^ 1 dp/drj. 
Taking these new variables and equations (|27[) and (|33[) . we 
find the dimensionless fe|| and fcj_ as, 



fcj 



fc||.R : 



k^R = 



a 2 p 



-ypC s cos 2 (q z ) 

r 2 n 1/2 
a p 



1 

m 2 



-,1/2 



b 2 



(37) 



(38) 



and we define also the dimensionless acoustic critical fre- 
quency, a ac , as 

Cd 



2H 



cos(a z ), 



(39) 



where Cd is the dimensionless sound speed defined as Cd = 

(c s7 p/p) 1/2 . 

With these new variables we derive expressions for the 
dimensionless slow and fast wave solutions and combine 
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them to obtain the expressions for the dimensionless paral- 
lel, en = £\\e at , and perpendicular, e± = e±e CTt , components 
of the displacement. For oscillations with frequencies larger 
than a ac , we get 



As il — k nri + o-t) ftAfbxbz 



ex = A; J (2Vw)e lCTi + 



/3A s <jp &a; 

7P 3/2 C S bl 



2U 



1 + 2iHfci 



H-kuri + at) 



(40) 



(41) 



and for oscillations with frequencies smaller than cr ac , we get 

(42) 



-As i(fc|,77+o-t) /3Afb x b z . — 



ex = A / J () (2 v /xp)e* ,Tt + 
PA s ap b x 2H 



7P 3/2 C S 6 Z 



i(fe|iij+CTt) 



(43) 



where A s is a complex, depth independent amplitude which 
depends on latitude, defined in terms of the previously 
adopted amplitudes by A s = A/(Ri/po) or A s — B/(Rs/po), 
depending on the colatitude. Similarly, Af is a complex 
depth independent amplitude which depends on colatitude, 
X = H 2 o 2 1 |6o| 2 , and t is the dimensionless time. 

The complex amplitudes A s and Af can be written 
in terms of their modulus, \A S \ and \Af\, and phases, <f) 3 
and <f)f. Before we proceed to calculate the disk integrated 
line-of-sight projected velocity, we still need to determine 
these amplitudes and phases. Their relative values are not 
arbitrary. Instead, they are imposed by the magnetoacous- 
tic coupling that takes place below, in the region where the 
magnetic and gas pressure are comparable. Thus, in order to 
find their values we fit our expressions for the displacement 
to th e numerical solutions that come out of the MAPPA 
code (|Cunhall200(l . From the fit, we extract the amplitudes 
and phases and, simultaneously, confirm that our analytical 
solutions for the displacement represent well the correspond- 
ing numerical solutions, in the limit when both f3 is much 
smaller than one, and the atmosphere becomes isothermal. 

Generally, since in the region of interest, /3 is much 
smaller than one, we may expect that the second term on 
the right hand side of equations (|40p - (|43|) can be neglected. 
That assumption, however, can break down if the relative 
value of the amplitudes \A S \ and \Af\ is such as to compen- 
sate for the small value of ft. We checked whether that was 
the case and found that the only situation in which the two 
terms become of the same order of magnitude is in equa- 
tion (|40|) . for colatitudes very close to 90°. In those cases, 
we have included the second term in our computations. In 
all other cases the second terms on the right hand side of 
equations (|40 [l - (|43p were neglected. 



4 THEORETICAL RADIAL VELOCITY 

In the observations what is seen is an integral over the visi- 
ble stellar disk. Therefore, we shall integrate the line-of-sight 
projection of the velocity over the visible stellar disk taking 




Figure 2. Schematic representation of the star. B is the local 
magnetic field vector, a is the angle between the magnetic field 
axis and the local magnetic field vector and S the angle between 
the magnetic field axis and the direction of the observer. 



into account the limb darkening effect, to simulate what ob- 
servers would see if they had the means to measure the radial 
velocity as function of atmospheric depth. 

In Fig[2]we show a schematic representation of the star, 
where the observer is considered to be at a general position. 
In the figure 8 is the angle between the magnetic field axis 
and the direction of the observer, which we shall assume to 
be fixed, meaning that we are considering the observer's view 
at a particular rotation phase. Moreover, B is the local mag- 
netic field vector and a is the angle between the magnetic 
field axis and the direction of the local magnetic field. In our 
analysis we consider two spherical coordinate systems: The 
spherical coordinate system, (r, 9, ip), aligned with the mag- 
netic field axis, and the spherical coordinate system aligned 
with the direction of the observer, (r, 8' ,ip'). The coordinates 
in the two sys tems are related through standard expressions 
given, e.g., by iKarttunen et al l l| 19961 ) . 

Taking the real part of equations (|40p and (I42f> . after 
neglecting the second terms on the right hand side, and dif- 
ferentiating with respect to time, we can write a general 
expression for the component of the velocity parallel to the 
direction of the magnetic field in the form, 



(-(Je|fc|| 1 17) 



sin(— qv|fc||j?7 + at), 



where q e and q r are factors defined in the following way: 
when the oscillation frequency is larger than the acoustic 
critical frequency, q e — and q r = 1; when the oscilla- 
tion frequency is smaller than the acoustic critical frequency, 
q e — 1 and q r — 0. Similarly, from equations (|4ip and Q43p . 
after neglecting the second terms on the right hand side, we 
find that the component of the velocity perpendicular to the 
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direction of the magnetic field is given by, 
v±(j),6,t) = -aAfJ (2^fxp) sin(fft). 

Assuming a linear limb darkening law, with a as the 
limb darkening coefficient, the expression for the velocity 
component parallel to the line of sight averaged over part of 
the visible stellar disk for a general position of the observer, 
Vint , thus becomes, 

sin(-gr|fe|||f? + <rt)^|| - A f J (2y/xp) sm(at)X±~\ 

x -^-(1 - o(l - cos0')) sin(26')d6'd(p', (44) 

where ip'f, 9[ and #^ are, respectively, the limit in lon- 
gitude and latitude, in the observer's spherical coordinate 
system, of the region of the visible stellar disk over which 
the integration is to be made. For an integration over the 
whole visible stellar disk 0' varies within the interval [0, n/2] 
and tp' varies within the interval [0, 27r]. Moreover, C n is a 
normalization factor defined as, 



ft rt 



(1 - o(l - cos 0')) sm(20')d0'd<p', 



'v'i Je ', 

and Xii and X± are the line-of-sight projections of the unit 
vectors e|| and e±, respectively. 



4.1 Fitting to Acos(ert+<?!>) 

In spectroscopic observations the radial velocity derived 
from the time series analysis is usually fitted to a function of 
the form Acos(crt+</>). From this fit, the autho rs derive the 
ampl i tude and the phase of the oscilla t ions (e.g. Kurtz et al.l 
120051 . l2006bl ; iRvabchikova et al.ll2007l ; ISachkov et al.ll2007l ). 
Since our ultimate goal is to compare the results of our work 
with the observations, in our study we fit the average line- 
of-sight velocity, Vint, to a function similar to the one above, 
and derive the amplitudes and phases of the oscillations as 
function of height in the atmosphere. From the fit we find 
the expression for the fitted phase, <f>, namely, 



<j>{S,ri) = arctan([ j** ^ ' [C M cos(-g r |fe|| \r, + ^e^ 1 '" 



+ C±cos(ct>f)]de'dip'] X [y J ' {-C u sm{-q r \k u \ri + <p B ) 

e -<ie\k u \v _ C± S i n (0 / )]de' V] _1 ) , (45) 
and the expression for the fitted amplitude, A, namely, 



v2C n 



[[ 



[— C| sin(— ?r|fc||h + <t>s) 



-9e|fe||l»J 



n 2 r r e 'f rv'f 

■ CjLBtn{<i> s )\do'd<p'\ +yj J [C 



-9el*||ll 



1-i 1/2 



cos(— 3 r |fe|||7y + <t> s ) + C± cos(0/)]d0'dy;'j 
where Cm and Cj_ are given by, 



(46) 



C|| = ^(l-a(l-cos0O)siii(20')*||: (47) 



and 



C*x = \A f \J (2^/xp)(l-a(l-cose'))sia(29')X±. (48) 



4.2 A toy Model 

In sec. [5] we will use expressions (|45[) to ((48} to study the 
behavior of the phase and amplitude of Vi n t in a number of 
case studies. Here we try to anticipate those results through 
an analytical analysis of the sum of wave-like functions. This 
simple analysis will help support the interpretation of the 
results obtained in the case studies. 

The motivation for the form of the wave-like functions 
that we shall consider comes from equation (|44[) . In the case 
of an observer aligned with the magnetic field axis, this equa- 
tion can be written as the sum of three integrals, namely, 



Vint (6, i], t) = h + Ii + h, 
where, 



h= 



-2- e" 1 k N 1 " sin{at +<f> s ) C\ \ dO'dip' , 

On 



7 a 



— s'm(—\k\\\r]-\-at- i r(t) s )C\\d8'd(p', 



2 C n 



i((Tt+ct>f)C±de'dip', 



(49) 



(50) 



(51) 



(52) 



where 8' c = ac is the critical angle above which the acoustic 
component will start to propagate as a running wave. Thus, 
I\ contains the information of the evanescent acoustic solu- 
tion only, 72 contains the information of the running acoustic 
solution only, and 73 contains the information of the stand- 
ing magnetic solution only. In the above we have assumed 
O'i ^ Q'c ^ Of. When 9' c < 8'i, Ii is zero and the lower colati- 
tude limit in 72 is taken to be 9[. Like wise, when 0' c > 0'f, 
I2 is zero and the upper colatitude limit in 7i is taken to be 



4-2.1 Visual superposition of waves of different nature 

In this toy study we define three simple waves, with veloc- 
ities vi, V2 and V3, inspired, respectively, by 7i, 72, and 73. 
The velocities are assumed to have the following form, with 
the subscripts r and i indicating the real and imaginary 
parts of the wavenumber, respectively, 



Vj = ipj sin(—kj r rj + at + <j>j), 



(53) 
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with j £ {1,2,3}, and amplitudes and wavenumbers char- 
acterized as follows: 



$1 =4, 1 e v/2H e- kl > rl ■ k u 



ki 



ki r = 0, 



f/>3 = ipse 



k 3i 



kii = ; k2 r — k2 



k 3r = 0, 



(54) 



(55) 
(56) 



where ip\, -02, tps, k\, k2 are all real, constant, and positive. 

Equations (|53[) and (|54[1 define a wave-like function that 
retains the characteristics of the integral of the line-of-sight 
acoustic component of the velocity when the oscillation fre- 
quency is smaller than the acoustic critical frequency, equa- 
tions (|53[) and ((55} define a wave-like function that retains 
the characteristics of the integral of the line-of-sight acous- 
tic component of the velocity when the oscillation frequency 
is larger than the acoustic critical frequency, and equations 
()53[) and (|56[) define a wave-like function that retains the 
characteristics of the integral of the line-of-sight magnetic 
component of the velocity. 

With these three simple wave-like functions in hand, 
we start by considering the simple case in which one of the 
contributions can be neglected when compared to the other 
two. In that case we define our total velocity as, 



Vtoy = Vj + V k 



(57) 



where j and k £ {1,2,3} are indexes that depend on the 
type of contribution that we are considering in the sum. 

Similar to what we did to Vint and to what is done 
in high-resolution spectroscopic observations, we fitted the 
equation for the velocity, v toy , to an expression of the type 
Atoy cos(at + (j)toy), and from the fit we derive the expression 
for the phase, 

<j>toy = arctan cos(—kj T -q + <j>j) + ip k cos(-k kr n + (f> k )) 

sin(-fc Jr ?7 + <pj) - Tp k sin(-fc fcr 77 + ^fc)) -1 ] , (58) 
and for the amplitude, 



Atoy = {ipj-ipk) +4V>jV>lb< 



-(k jr — k kr )ri + (f>j — 4>k 



.(59) 



These expressions are analogous to equations (I45p and 
(|46[) . respectively, except that they are substantially simpler 
because they do not contain the latitudinal dependence of 
each of the integral components. 

Concerning the amplitude, Atoy{if), we can see it will 
be affected by the relative size of the amplitudes tpj and 
$ k and by the phase [ ~ (fcj> d^H±hzl ^ j . if ^ i s mucn 

larger, or much smaller, than if> k , the first term on the right 
hand side of equation ()59|) will dominate and no structure in 
the amplitude as function of atmospheric height is expected. 
Thus the behaviour of the latter with atmospheric height 
will be similar to the behaviour of the predominant com- 
ponent. That is the most common situation. On the other 
hand, if the two components have similar amplitudes, i.e., 
ipj ~ ipk, then the second term on the right hand side of 
equation (|59p may dominate and some structure in the am- 
plitude as function of atmospheric height may be expected. 
However that structure will depend strongly on the phase, 
— kk rJi+4> J — tfej ^ j n region under consideration. In 



fact, when values typical of those found in our general anal- 
ysis are taken for the latter, we find that in most cases no 
structure at all is found, i.e., the amplitude simply increases 
with atmospheric height as in the first situation considered. 

The results for the phase, (ptoy, will also depend on the 
relative amplitude of each component and on the argument 
of the trigonometric functions. However, in this case no mat- 
ter how large the value for the amplitude of one of the com- 
ponents is, if its phase is such that either the cosine or the 
sine of the corresponding argument goes through a zero, the 
component with smaller amplitude will be the one dominat- 
ing the corresponding sum (nominator or denominator) in 
equation (|58l) . This, as we will see later, will result in a di- 
versity of phase behaviour that will be possible even when 
one component is much larger than the others. 

For completeness, and because this analytical analysis 
will be of relevance for the interpretation of the results pre- 
sented in section below we provide the expressions for the 
phase and amplitude obtained in the case of adding three 
simple wave-like functions, v toy = Vj 
we find, 



'Jtoy 



= arctan 



^ipj cos(-k Jr f] + 4>j) + ip k cos(-k kr n + tf> k ) 



+ip p cos(—k Pr ri + (j) p )j ^ — ipj sm(—kj r r] + 4>j) 
-4>k sin(-fcfe r ?7 + <f>k) + tjjp sm(- k Pr rj + <j> p )) j, (60) 



and 



A 2 - 



(i>j sxa(-kj T rj + 4>j) + tjj k sin(-k kr r] + 4> k ) 
i> p sin(— k Pr rj + (p p )J + (jjj cos(— kj r r) + 4>j) 



+ 



+ 



ijj k cos(-fc fer ?7 + <f> k ) + 4> p cos(- k Pr r] + <j> p ) ) . (61) 



4-2.2 Visual superposition of waves of similar nature 



Based on equation (|49|) we have suggested three components 
for the velocity, inspired, in particular, in the two solutions 
for the acoustic component of the velocity and in the so- 
lution for the magnetic component. However, if we look at 
the arguments of expressions (I50[) . (|51[) and (|52[1 . we no- 
tice that the amplitudes, the vertical wavenumbers and the 
phases all depend on colatitude. This was not considered in 
section 14.2.11 because our intention was to eliminate the in- 
tegrals, to keep the analysis simple. Nevertheless, it is clear 
that such dependence can itself introduce a variety of phase 
behaviour, even when a component of a particular nature is 
considered alone. To explore this fact, we can associate Vj 
and v k in equation ((57} to two components of similar nature, 
with different amplitudes, vertical wavenumbers and input 
phases. 

As an example, let us consider a single evanescent wave. 
The amplitude of that wave increases with atmospheric 
height, while the phase stays constant. On the other hand, if 
we have a superposition of two evanescent waves with differ- 
ent wave numbers and phases, we can find a phase varying 
with atmospheric height. Similarly, if we consider a single 
outwardly running wave, the amplitude increases and the 
phase decreases with height. However, if we have a superpo- 
sition between outwardly running waves with different wave 
numbers and phases, we can have a phase scenario different 
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B p (kG) 


(mHz) 


Limits of integration 


Case 1 


4 


1.0 


whole stellar disk 


Case 2 


4 


2.5 


whole stellar disk 


Case 3 


1 


1.7 


whole stellar disk 


Case 4 


1 


2.1 


whole stellar disk 


Case 5 


2 


3.0 


whole stellar disk 


Case 6 


4 


2.2 


[62°,85°] 



Table 1. Characteristics of the Case Studies. All cases are for an 
observer pole-on and 1 = 1. 



from the latter. Finally, when considering a single standing 
solution defined by equations (|53(l and (I56|l , both amplitude 
and phase are constant with atmospheric height. In this case, 
if we have a superposition of two of these solutions with dif- 
ferent amplitudes and phases, the only thing that may be 
affected is the sign of the phase, since in all cases the total 
amplitude and phase will remain constant. 



5 RESULTS AND ANALYSIS 

In this section we present a selection of the results obtained 
in our study of the integrated line-of-sight velocity. The 
cases to be presented were chosen so as to illustrate dif- 
ferent behaviours of the phase. However, as our aim was 
to conduct an in-depth study of these cases, we have re- 
stricted the parameter space considerably. In particular, we 
have fixed the position of the observer, which is always con- 
sidered to be aligned with the magnetic axis, and we have 
considered modes which in the absence of a magnetic field 
would be pure dipole modes (I = 1). The results of a system- 
atic study of the parameter space, including variations of 5 
and I, conducted with the same approach, will be presented 
in a separate paper. A list of the case studies is shown in 
Table [fl 

In the computations we have used a stellar model taken 
from the Code d'Evol ution Stellaire Adaptati f et Modulairc 
(CESAM code), (e.g. iMorel fc Lebretorj|20og) . with a mass 
M = 1.8 Mo, a radius R = 1.57 R©, and an effective tem- 
perature T e g = 8363 K. To this model, we have matched 
an isothermal atmosphere at a temperature Ti so = 6822 K. 
The acoustic cut-off frequency of the model is 2.7 mHz. We 
have considered different values for the polar magnetic field, 
Bp, different cyclic oscillation frequencies, u, varying from 
1.0 mHz to 3.0 mHz, and different integration intervals, to 
take into account the possibility that the lines from which 
the radial velocities are derived in the observations are not 
homogeneously distributed over the stellar surface but, in- 
stead, concentrated in particular regions of the star. 

For each case we derived the amplitude, A, and the 
phase, (j>, of the integrated line-of-sight velocity, using equa- 
tions l|45p and ([46]) , respectively. Moreover, we fitted equa- 
tion (|50|) to a function of the type A ev cos(at + 4>ev) to derive 
the amplitude, A ev , and the phase, (j> ev of the evanescent 
component; we fitted equation (|51|) to a function of the type 
A r cos(crt + 4>r) to derive the amplitude, A r , and the phase, 
4> r of the running component; and we fitted equation (|52p to 
a function of the type A„t cos (erf + (j> s t) to derive the ampli- 
tude, A st , and the phase, 4>st of the standing component. In 





Amplitude 


Phase 


Nature 


Total velocity 


A 


4> 


Ac + Mag 


Evanescent component 




<j>ev 


Ac 


Running component 




4>r 


Ac 


Standing component 


A st 


4>st 


Mag 


Total acoustic component 


A ac 


4>ac 


Ac 



Table 2. Different representations for the amplitude (second col- 
umn), phase (third column), and nature of the wave (where Ac 
stands for acoustic and Mag stands for magnetic), according to 
the component of the velocity from which they were derived (first 
column). 



addition, we fitted the total acoustic term of Vi n t to a func- 
tion of the type A ac cos(<rf+</> ac ) to derive the corresponding 
amplitude, A ac , and phase, 4> a c- In Table [2] we show the dif- 
ferent representations of the amplitude and phase according 
to the component of the velocity integral from which they 
were derived. 

In the plots below, the amplitudes and phases are shown 
as a function of atmospheric height, given in units of the 
dimensionless scale height. We note that the zero of the at- 
mospheric height corresponds to the base of the isothermal 
atmosphere. In practice, the region where rare earth lines 
are formed in roAp stars will be only a small slice of the 
layers that we are considering. However, the phase of each 
component contributing to iw in that region will depend 
on the model considered, as well as on the magnetic field 
intensity. Since exploring a large grid of magnetic models is 
out of the scope of this work, we extend the isothermal part 
of our model atmosphere for 10 scale heights, in the hope of 
finding a variety of combinations of 4>st, <f>r, an d <f) ev , and, 
thus, all possible behaviours of the total phase. 

5.1 Pulsation Amplitude 

Concerning the amplitude, we found that in most cases 
the latter increases with atmospheric height. In Fig. [3] we 
show an example of this behavior, for the case of B p =2 
kG, v = 3.0 mHz, and an integration over the whole vis- 
ible stellar disk. This behavior for the amplitude was to be 
expected given our findings in the toy study presented in 
section 14.21 and can be explained by the fact that the den- 
sity is decreasing with height. In this case, the amplitude 
goes with 1/y/p. This is consistent with the running acous- 
tic component dominating the integral of the velocity. In 
Fig. [4] we show the amplitude derived for the case of B p =4 
kG, v = 1.0 mHz, and an integration over the whole visible 
stellar disk. While the amplitude in this case also increases 
with atmospheric height, it does so at a smaller rate than 
in the previous case. That is a consequence of the integral 
for the velocity being dominated by the evanescent, rather 
than the running, acoustic component. 

The only exception to this growing behaviour of the 
amplitude that shall be shown in the present paper is that 
discussed in the case study 6, presented below. Other ex- 
ceptions are found when the parameter space is not as re- 
strictive as the one considered here. Those cases shall be 
discussed in the upcoming paper. 
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atmospheric height 

Figure 3. Dimensionless amplitude as a function of atmospheric 
height (in units of H) - B p =2 kG, v = 3.0 mHz, and for an 
integration over the visible stellar disk. 
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Figure 4. Dimensionless amplitude as a function of atmospheric 
height (in units of H) - B P =A kG, v = 1.0 mHz, and for an 
integration over the visible stellar disk. 



5.2 Pulsation Phase 

Concerning the phases, we find a situation that is much more 
diverse. This was also to be expected given our findings in 
section 14.21 

Case 1: Constant phase 



In case 1 we consider a magnetic field with polar mag- 
nitude B p =l kG, and a cyclic oscillation frequency v — 
1.0 mHz. Moreover, the integration is over the whole visible 
stellar disk. The oscillation frequency in this case is smaller 
than the acoustic cutoff frequency. The critical angle, 9 ac , 
is 78°, which means that for colatitudes smaller than this 
angle the acoustic term will be an evanescent wave, with a 
phase independent of atmospheric height, while for larger 
colatitudes the acoustic term will be an outwardly running 
wave with phase decreasing with atmospheric height. 

As we can see in Fig[5] the phase obtained with equation 
(|45p is constant throughout the region of the atmosphere un- 
der consideration. In Fig. [6]we show the amplitudes derived 
from each of the components of the velocity. Clearly, the 
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Figure 5. Case 1 - B p =4 kG, u = 1.0 mHz, and integration 
over the whole visible stellar disk. Total phase overlapped with 
phase derived from the evanescent term alone as a function of 
atmospheric height (in units of H). The two lines are coincident. 
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Figure 6. Case 1 - same as in figure[5] Left panel: Dimensionless 
amplitudes derived from each term of the velocity: full line - A ev ; 
dashed line - A r ; dashed-dotted line - A B t Right panel: Zoom in 
the y axis for A r . 



amplitude derived from the evanescent term alone is always 
much larger than the amplitudes derived from the running 
or magnetic terms alone. 

In Fig[5] we also show the phase derived from the 
evanescent term alone. It overlaps with the line for the total 
phase. This is a simple case, the simplest of all cases that 
we are going to consider here, in which one of the terms 
of the velocity dominates entirely the final result, due to its 
large amplitude when compared to the amplitudes exhibited 
by the other terms. We note, however, that having a dom- 
inating amplitude is not a sufficient condition to dominate 
entirely the final results, as we will see in other examples. 

It is important to notice that other situations exist in 
which the resulting phase is constant due to reasons different 
to the one considered in this case study. If evanescent or 
standing waves dominate the integral of the velocity, and 
if the phases derived from the corresponding terms alone 
are constant, the resulting phase will generally be constant. 
Moreover, the resulting phase can sometimes be constant 
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tion over the visible stellar disk. Total phase (full line) and phase 
derived from the acoustic term alone (dashcd-dottcd line) as a 
function of atmospheric height (in units of H). 



if the visual superposition happens between evanescent and 
standing waves with similar amplitudes. 

Case 2: Phase increasing deeper and decreasing higher 
in the atmosphere 

The case that we present here is for Bp—4 kG, v = 2.5 
mHz, and an integration over the whole visible stellar disk. 

This case is similar to the previous one in the sense that 
the oscillation frequency is smaller than the acoustic cutoff 
frequency. However, here the critical angle is only 43°. The 
phase derived for this case is shown in Fig. (full line). 
Since the phase behaves differently, in different regions of 
the atmospheric layers, we will analyze the outermost layers 
and innermost layers separately. 

Higher in the atmosphere, in the region where the phase 
decreases with atmospheric height, the amplitude derived 
from the evanescent term alone is again much larger than 
the amplitude derived from the running or magnetic terms 
alone, as can be seen in Fig. [8] 

Thus, in the outermost layers this case seems similar 
to case study 1, the main difference being that in the lat- 
ter the total phase was constant, while here it decreases 
slightly with atmospheric height. To understand why the 
phase varies with atmospheric height, we will first consider 
the phase of the dominating term, namely, cj> e v, which is 
shown on Fig. [§] (full line). 

Unlike what was found in case study 1, in the present 
case <j>ev varies with atmospheric height. To check that we 
understand the reason for this behaviour, we turn to the toy 
study presented in section 14.21 In our toy study the anal- 
ogy with this case is achieved by considering two evanescent 
waves. The amplitude of the first wave, tpi, is considered to 
be larger than the amplitude of the second wave, ip2, inspired 
in the fact that the amplitude of the acoustic component of 
the velocity decreases with colatitude. Moreover, the acous- 
tic wave number of the first wave, ki, is considered to be 
larger than the acoustic wave number of the second wave, 
/C2, inspired in the fact that the absolute value of the acous- 
tic vertical wavenumber for the evanescent term decreases 
with colatitude. 
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Figure 8. Case 2 - same as in figurefT] Left panel: Dimensionless 
amplitudes derived from each term of the velocity: full line - A ev ; 
dashed line - A r ; dashed-dotted line - A s t Right panel: Zoom in 
the y axis for A r . 
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Figure 9. Case 2 - same as in figure Phase derived from the 
evanescent term alone (full line) and from the running term alone 
(dashed-dotted line) as a function of atmospheric height (in units 
of H). 



The contributions of the two waves to the numerator 
and denominator of equation (|58|) depend critically on the 
phases (j>i an d cf>2- These can be estimated from the input 
phases of our real model. In Fig. [TU]we show the input phase 
of the acoustic component, ci s , in the region of colatitude 
where the acoustic component is evanescent. We conclude 
from the plot that the value of <f>i (the wave nearest the 
magnetic pole) is such that | cos(<^>i)| is close to 1 and thus 
larger than | cos(02)|- Since the amplitude, ipi, is also larger 
than ip2 we are going to assume that the contribution of 
the second wave to the numerator of equation (|58[) can be 
neglected, in which case the latter becomes, 



arccot(— tan(0i) 



l\>2 



(ki-k 2 )?7 ' 



(62) 



The behaviour of 4>toy will therefore depend critically 
on the sign of cos(oii) and sin(cS2). From Fig. [10] we see 
that cos(<f) s ) is always negative and, thus, so is cos(cii), and 
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Figure 10. Case 2 - same as in figure [7] Zoom of the input 
phase of the acoustic term, <p s , in the region where the solution 
is evanescent as function of colatitudc. 



< 
< 



0.0010 




0.0005 



0.0000 

0.0 0.5 1.0 1.5 2.0 

otmospheric height 
Figure 11. Case 2 - same as in figure [7] Zoom in the xx axis 
for the dimensionless amplitudes derived from each term of the 
velocity: full line - A e „; dashed line - A r ; dashed-dotted line - 
A s t. 



sin((j} s ) is always positive and, thus, so is sin(<jf>2). There- 
fore, since e^ kl ~ k2 ' lv increases with atmospheric height, 4>t oy 
will decrease with atmospheric height. We note that such 
conclusion would be reached if the term of </>2 in the de- 
nominator had not been neglected, but the algebra would 
be more complicated. 

Considering the results of our model, and equation (|62[) . 
we may be tempted to conclude that the phase behavior 
seen in the outermost layers in this case study is due to 
a visual superposition of evanescent waves of different wave 
numbers and input phases, with the resulting phase decreas- 
ing with atmospheric height. However, when comparing the 
phase derived from the evanescent term alone, (full line in 
Fig©, with the total phase shown in Fig.[7]higher in the at- 
mosphere, we find that their slopes are not the same. There- 
fore, one of the other two components of the velocity must 
be contributing to the behavior of the total phase, fn fact, 
looking at Fig. [7] where the phase derived from the acoustic 
term alone is shown (dashed-dotted line), one can see that 
the slope of the latter is very similar to the slope of the 
total phase in the outermost layers. Thus we may conclude 
that it is the acoustic running component, whose phase also 
decreases with atmospheric height (cf. Fig©, that is influ- 
encing, together with the evanescent component, the behav- 
ior of the phase in the outermost layers of the atmosphere. 
This is an example of a situation in which a component of 
significant lower amplitude influences the final result. 

Deeper in the atmosphere, in the region where the phase 
is increasing with atmospheric height, the amplitude derived 
from the evanescent term alone is still larger than the am- 
plitude derived from the running and magnetic terms alone, 
as can be seen in Fig. 1111 Thus, we could naively expect the 
evanescent term to dominate the phase behavior. However, 
from a comparison of Figs[7]and[9l we see that such is clearly 
not the case. 

Looking back at our toy study, the analogy with the 
case presented here is achieved by considering the first wave 
to be an evanescent wave, the second wave to be a run- 
ning wave, and the third wave to be a standing wave with 
constant amplitude. The amplitudes of the wave-like func- 
tions in the toy model are estimated from the amplitudes 
derived when considering the evanescent, the running, and 
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Figure 12. Case 2 - same as in figure[7] Phase derived from the 
magnetic term alone as a function of atmospheric height (in units 
of H). 



the standing terms alone. Similarly, the phases in the toy 
model are estimated from the phases of the corresponding 
components, adding tt/2, i.e., <j>i = e „ + 7r/2, (j>2 = </> r + 7i"/2, 
and 03 = <f>st + 7r/2. The 7r/2 difference is a consequence of 
using the sin function in equation (|53[1 and the cos function 
in the fits to extract <f) ev , <f) r , and <f) s t. 

Due to its amplitude, we find, after checking the terms 
in equation (|60[1 , that the contribution of the running wave 
to (f> t0 y can be neglected. Moreover, from the inspection of 
(f> s t we find that we can also neglect the contribution of the 
standing wave to the denominator of equation (|60[) . How- 
ever, the contribution of the latter to the numerator of equa- 
tion (|60[l cannot be neglected, as happened higher in the 
atmosphere. In fact, the amplitude of the evanescent wave 
deeper in the atmosphere is much smaller than its ampli- 
tude higher in the atmosphere, and when the phases of the 
evanescent and standing components are taken into account, 
the corresponding terms in the numerator of equation (|60|l 
become of the same order of magnitude. Then equation (|60|l 
becomes, 
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Figure 13. Case 3 - B p =l kG, v = 1.7 mHz, and an integration 
over the visible stellar disk. Total phase as a function of atmo- 
spheric hcight(in units of H). 



, . , «. i>3 v/2-H -k lv COs((j>3) 

(fitoy = arctanl — cot(ffli) — e" e . . , , ). (63) 

The behaviour of (fitoy will, therefore, depend critically 
on the sign of cos(03) and sin(<jf>i). From Fig |12l and the rela- 
tion 03 = (f) s t +7r/2, one concludes that cos(03) is negative, 
and from Fig. [9] and the relation cfii = cfi av + n/2, one con- 
cludes that sin(0i) is positive. Since, e v ^ 2n e~ kiv increases 
with atmospheric height, and <fi\ varies only very little within 
those lower layers, the result will be a phase increasing with 
atmospheric height. 

Considering the results of this case study, and equa- 
tion (|63p . we may thus come to the conclusion that the phase 
behavior seen in the deeper layers is likely due to a visual 
superposition of the evanescent acoustic waves and magnetic 
waves, with the resulting phase increasing with atmospheric 
height. 

Case 3: Phase increasing with atmospheric height 



Case study 3 is for B p =l kG, v — 1.7 mHz, and an 
integration over the whole visible stellar disk. The oscillation 
frequency in this case is smaller than the acoustic cutoff 
frequency. The critical angle is 67°. 

From Fig |13l we see that in this case the phase obtained 
with equation (|45[) is increasing throughout the region of 
the atmosphere under consideration. In Fig. [TJwe show the 
amplitude derived from each of the components of the ve- 
locity alone. As before, we find that the amplitude derived 
from the evanescent term alone is always much larger than 
the amplitude derived from the running or magnetic terms 
alone. Yet, the behaviour of the phase is different from that 
of the case studies considered before. 

To understand the phase behaviour, we again look for 
the analogous case in our toy model, taking the first, second, 
and third waves, as evanescent, running, and standing waves, 
respectively. 

Due to its amplitude, we find after checking the magni- 
tude of the terms in equation (|60|) . that the contribution of 
the running wave to (fitoy can again be neglected. Moreover, 
looking at Fig |15l where we show the phases derived from 
the standing and evanescent terms alone, respectively, and 
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Figure 14. Case 3 - same as in figure [131 Left panel: Dimension- 
less amplitudes derived from each term of the velocity: full line 
- A e i>; dashed line - A r ; dashed-dotted line - A s f Right panel: 
Zoom in the y axis for A,-. 
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Figure 15. Case 3 - same as in figure [T31 Phase derived from the 
evanescent term alone (full line) and phase derived from the mag- 
netic term alone (dashed-dotted line) as a function of atmospheric 
height (in units of H). 



taking into account that 03 = <fi st +ir/2 and <f>\ — </> ev +Tv/2, 
we see that we can neglect the contribution of the standing 
wave to the numerator of equation (|60[) . However, the con- 
tribution of the latter to the denominator of equation (|60|) 
cannot be neglected, since sin(0i) is sufficiently small that 
the evanescent and standing contributions become of the 
same order. Then equation (|60p becomes, 

i , % ^3 -, ; /2« k 17 ,Sm(03)s . . 

(fitoy = arccot(-tan(0i) — e " e 1 ' 7—7)- (64) 

-01 cos(0i) ; 

From Fig. 1151 one infers that sin(03) is positive and 
cos(0i) is negative. Since, e~ v ^ 2H e hlv decreases with atmo- 
spheric height and tan(0i) varies only very little within the 
atmospheric region considered, the result will be a phase 
increasing with atmospheric height. 

From equation (|64[) we may thus conclude that the 
phase behavior seen in this case study is likely due to a 
visual superposition of the evanescent and magnetic terms, 
with the resulting phase increasing with atmospheric height. 
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Figure 16. Case 4 - B p =l kG, v = 2.1 mHz, and an integration 
over the visible stellar disk. Total phase (full line), phase derived 
from the evanescent term alone (dashed-dotted line) and phase 
derived from the magnetic term alone (dashed line) as function 
of atmospheric height (in units of H). 
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Figure 17. Case 4 - same as in figure [16] Left panel: Dimension- 
less amplitudes derived from each term of the velocity: full line 
- A ev ; dashed line - A r ; dashed-dotted line - A s t Right panel: 
Zoom in the y axis for A r . 



We note that there are other situations that can result 
in the phase increasing with atmospheric height. Such be- 
havior could be found as a result of a visual superposition of 
running and evanescent waves, or of running and standing 
waves, under particular conditions. 

Case 4: Phase decreasing deeper and constant higher 
in the atmosphere 



The present case is for B p — 1 kG and v = 2.1 mHz. The 
integration is over the whole visible stellar disk. The oscilla- 
tion frequency is smaller than the acoustic cutoff frequency 
and the critical angle is 57°. As in case study 2, in this case 
the phase behaves differently, in different atmospheric layers. 
Thus, we will analyze the outermost layers and innermost 
layers separately. 

Higher in the atmosphere, in the region where the phase 
is approximately independent of the atmospheric height, the 



amplitude derived from the evanescent term alone is again 
significantly larger than the amplitude derived from the run- 
ning or magnetic terms alone, as we can see in Fig |17l More- 
over, if we compare the evanescent phase shown in Fig |16l 
(dashed-dotted line), with the total phase in the same fig- 
ure (full line), we may conclude that the result is consistent 
with the evanescent term dominating the phase behaviour 
in that region of the atmosphere. The difference between the 
values of the phase can be explained by the influence of the 
lower atmospheric layers on <f>. 

Deeper in the atmosphere, in the region where the phase 
is decreasing with atmospheric height, the amplitude derived 
from the evanescent term alone is still larger than the ampli- 
tudes derived from the running and magnetic terms alone. 
However, looking at the phase derived from the evanescent 
term alone, we can see that it does not show the same be- 
havior, deeper in the atmosphere, as the total phase. Thus, 
we may conclude that at least one of the other terms must 
be influencing the behavior of the total phase. 

To investigate further, we look back at our toy study, 
identifying, as before, the first, second, and third waves with, 
evanescent, running, and standing waves, respectively. Due 
to its amplitude and phase, the contribution of the running 
wave to 4>toy can be neglected. Moreover, from Fig. 1161 we 
see that cos(^)i) ~ cos{4> ev + tt/2) will not be small and, 
thus, we may neglect the contribution of the standing wave 
to the numerator of equation (|60|) . However, the contribu- 
tion of the latter to the denominator of equation (|60[) can- 
not be neglected, as happened higher in the atmosphere. 
As seen in Fig. 1171 the amplitude of the evanescent wave 
deeper in the atmosphere is considerably smaller than its 
amplitude higher in the atmosphere, and when the phases 
of the evanescent and standing components are taken into 
account the terms in the denominator of equation (160[) be- 
come of the same order of magnitude. Then equation (|60l) 
takes the form, 

4> toy = arccot(-tan(^) - ^ e -^* e ^f^M). (65) 

tpi COS (01 ) 

This equation is the same as equation (|64[) . However, in 
this case we deduce from Fig |16l that sin((j!>3) is negative and 
from Fig |16l that cos(^i) is also negative. Thus, in this case 
we find a phase decreasing with atmospheric height. 

From equation (|65|l . we may thus conclude that the 
phase behavior seen in the deeper layers in this case study 
is likely due also to a visual superposition of the evanescent 
and magnetic terms, but in this case the resulting phase 
decreases with atmospheric height. 

Case 5: Phase decreasing fast throughout the whole 
atmospheric layer 



This case study is for B p = 2 kG, v = 3.0 mHz, and an 
integration over the whole visible stellar disk. The oscilla- 
tion frequency is larger than the acoustic cutoff frequency, 
therefore we do not have a critical angle: the acoustic term 
will be a running wave at all latitudes. 

From Fig. 1181 we see that in this case the phase obtained 
with equation (|45p is decreasing throughout the whole region 
of the atmosphere under consideration. In Fig. [19] we show 
the amplitude derived from each of the components of the 
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Figure 18. Case 5 - B p =2 kG, v = 3.0 mHz, and an integration 
over the visible stellar disk. Total phase overlaped with phase 
derived from the running term alone as a function of atmospheric 
height (in units of H). The two lines are coincident. 
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Figure 20. Case 6 - B P =A kG, v = 2.2 mHz, and an integration 
of [62,85]° of colatitude. Dimensionless amplitude as a function 
of atmospheric height (in units of H). 
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Figure 19. Case 5 - same as in figure H"8l Dimensionless ampli- 
tudes derived from each term of the velocity: dashed line - A r ; 
dashed-dotted line - A.* 



velocity alone. Clearly, the amplitude derived from the run- 
ning term alone is always much larger than the amplitude 
derived from the magnetic term alone. 

Moreover, Fig. 1181 shows that the phase derived from 
the running term alone overlaps with the total phase, and 
we may conclude that the running term dominates also the 
phase behavior. This is a again a simple case, in which one of 
the terms of the velocity dominates entirely the final result. 

Nevertheless, it is important to notice, as before, that 
other situations can exist in which the resulting phase de- 
creases with atmospheric height due to reasons different 
than the one considered in this case study. Whatever the 
visual superposition, if it involves running waves that, due 
to their amplitude and phase, are dominating, the result 
will be a phase decreasing with atmospheric height. A phase 
decreasing with atmospheric height can also occur due to 
a visual superposition of running and evanescent waves, of 
running and standing waves, or even of evanescent waves 
with slightly different wave numbers. Despite this, we note 
that when the resulting phase is decreasing with atmospheric 
height due to a visual superposition of waves of different na- 
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Figure 21. Case 6 - same as in figure [201 Phase as a function of 
atmospheric height (in units of H). 



ture, the rate at which it decreases is significantly smaller 
than when it is due to the running waves alone. The only 
exception to this is the case when the pulsation amplitude 
approaches a zero, and the phase decreases significantly, in 
some limited region of the atmosphere. That is the situation 
that will be investigated in case study 6, below. 

Case 6: A node in the amplitude and a phase jump 

The last case study to be considered is for B p =4 kG, 
v = 2.2 mHz, and an integration over colatitudes within 
the interval [62°, 85°]. The oscillation frequency is smaller 
than the acoustic cutoff frequency and the critical angle is 
55.5°. Since the integration is over a colatitude region that 
is above the critical angle, the acoustic term has always the 
form of a running wave, with a characteristic wavelength 
and amplitude which depend on co-latitude. 

In Figs. 1201 and |2"T1 where the results for this case study 
are shown, we can see a zero in the amplitude accompa- 
nied by a phase jump of ~ it, at an atmospheric height 
close to AH. Although at first sight this behaviour may be 
thought to result from a node in the velocity eigenfunction, 
we know that such cannot be the case. In fact, from our 
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Figure 22. Case 6 - same as in figure 1201 Dimensionless acoustic 
component of the integrated line-of-sight velocity as a function of 
atmospheric height (in units of H). Different lines are for different 
instants during the pulsation cycle. 




atmospheric height 
Figure 24. Case 6 - same as in figure 1201 Dimensionless inte- 
grated line-of-sight velocity as a function of atmospheric height 
(in units of H). Different lines are for different instants during the 
pulsation cycle. 
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Figure 23. Case 6 - same as in figure[20] Dimensionless magnetic 
component of the integrated line-of-sight velocity as a function of 
atmospheric height (in units of H). Different lines are for different 
instants during the pulsation cycle. 
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Figure 25. Case 6 - same as in figure 1201 Dimensionless ampli- 
tudes derived from each term of the velocity: dashed line - A ac ; 
dashed-dotted line - A s t- 



study of the velocity field we know that in this region of 
our model atmosphere the component of the velocity along 
the magnetic field is, at each latitude, a running wave, and 
the component perpendicular to the magnetic field has an 
amplitude very closely independent of atmospheric height. 
Thus, the velocity eigenfunction cannot have the properties 
of a standing wave with a node in the atmospheric layers 
considered. We can, therefore, anticipate that the behavior 
seen in the amplitude and phase in this case study is only 
that of a "false node", that results from the visual super- 
position of the line-of-sight projections of the acoustic and 
magnetic components, when integrated over the stellar re- 
gion considered. 

In Figs. [22] and [23] we show, respectively, the acoustic, 
Vac, and magnetic, Vmag, components of the integrated line- 
of-sight velocity, v int . We note that in the present case v mag 
is identical to the integral ^3, defined in equation (|52[) . while 
Vac is a subsection of the integral I2, corresponding to the 
interval of colatitudes under consideration. Different lines 
correspond to different times, during a pulsation cycle. It is 
clear that none of the integrated velocity components shows 



the characteristics of a standing wave with a node. Never- 
theless, the signature of the apparent node is evident in Vint, 
shown in Fig. 1241 From the inspection of Fig. 1251 where we 
show the amplitudes derived from each of the integrated ve- 
locity components, we see that the apparent node results 
from the cancellation between the two terms of the integral 
for the total line-of-sight integrated velocity, when averaged 
over the region of the stellar disk considered. The jump in 
the total phase results from the change in the sign of Vi n t, 
like it would occur when crossing a node in a pure standing 
wave. 



6 DISCUSSION 

The results presented in the previous section are to be con- 
sidered as illustrative of what may take place in the outer 
layers of the atmosphere, hence, in the region where spec- 
tral lines associated with rare earth elements are formed. 
In deeper layers of the atmosphere, the assumed isother- 
mal condition becomes inappropriate and, depending on the 
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magnetic field intensity, the magnetoacoustic wave may not 
yet be decoupled into acoustic and magnetic components. 

If rare earths are indeed formed in the region of the at- 
mosphere considered in this work, then a phase jump of n 
associated with a zero of the pulsation amplitude, derived 
from these lines, cannot be due to a real pulsation node. 
Nevertheless, the same behaviour of the phase and ampli- 
tude derived from lines that are sensitive to deeper regions of 
the atmosphere (say, the base of the photosphere) may still 
reflect the presence of a real pulsation node. The analysis 
of the theoretical radial velocity near the base of the pho- 
tosphere is beyond the scope of this paper, thus we cannot 
make conclusions with regards to alleged nodes derived from, 
e.g, the in spection of the co re of H a lines, or of the Fe and 
Mg lines (|Elkin et alj|2005j ). However, there are reasons to 
believe that real nodes may exist at the phot ospheric level, at 
least in some latitudinal regions of the star l|Gautschv et al.l 
1 19981 ; iKhomenko fe Kochukhovll2009l ; ISaio et alj|20ich . and 
result in radial velocity amplitude and phase behaviours that 
are similar to those of false nodes. 

Despite the diversity of phase behaviour found in the 
results presented in sec. [5] it is clear that in the generality 
of the case studies considered the derived phase variations 
are rather small. In fact, some of these are unlikely to be 
detected, given the uncertainties in the phases derived from 
the time-series of spectra. Large and rapid phase variations 
are derived in case studies 5 and 6. The first of these concerns 
a situation in which the running wave largely dominates the 
integral for the radial velocity, while the second is the case 
of the fake node where the phase varies at a high rate in 
a localized region of the atmosphere. Additional situations 
of strong variations of the phase, as well as of false nodes, 
are found when different positions of the observer and mode 
degrees are considered. Those shall be detailed in paper II. 

Looking more generally at our results, we may conclude 
that the behaviour of the theoretically derived amplitudes 
is in good qualitative agreement with that inferred from the 
observations, the most common trend being an increase with 
atmospheric height. 

Concerning the pulsation phases, our findings illustrate 
well the difficulties that may arise when attempting to infer 
information about the properties of the pulsations in the 
atmospheres of roAp stars through the inspection of the 
spectroscopic time-series. Those difficulties are due, in par- 
ticular, to the observer's inability to resolve the star. The 
observations need to be thought of as visual superposition 
of different wave functions, and / or similar wave functions 
with different characteristic wavelengths, amplitudes, and 
phases, which are not always easy to disentangle. A clear 
example of this is the fact that the same general phase be- 
haviour can be explained by combinations of different wave- 
type solutions. 

In fact, phases that are independent of atmospheric 
height may indicate that evanescent waves of acoustic na- 
ture, or standing waves of magnetic nature, dominate the 
integral for the radial velocity, or that both type of waves 
contribute similarly to the latter. Similarly, different situa- 
tions can lead to an increase of the phase with increasing at- 
mospheric height, such as the visual superposition of evanes- 
cent and standing waves, the visual superposition of running 
and evanescent waves, or the visual superposition of running 
and standing waves. Moreover, this behaviour of the phase 



can also result simply from a visual superposition of evanes- 
cent waves with slightly different wave numbers. We should 
note, however, that for lines formed in the region of the at- 
mosphere considered in the present work, this behaviour of 
the phase cannot be due to a downwardly propagating wave, 
since no such type of solution is present there. Finally, phases 
decreasing with atmospheric height can result from running 
waves propagating upwardly, from a visual superposition of 
running and evanescent waves, or of running and standing 
waves, and from a visual superposition of evanescent waves 
with slightly different wave numbers. 

Despite the potential difficulties that may arise when 
attempting to disentangle the different possible scenarios, it 
is also clear that some of these lead to substantially different 
rates of phase changes, as well as amplitude changes with 
atmospheric height. Thus, it is likely that at least some of 
the apparent degeneracy in phase behaviour may disappear 
if a more realistic model of the stellar atmosphere, than 
that used in the present analysis, is considered, and a direct 
confrontation of observations and theoretical predictions is 
carried out. 

Finally, as acknowledged before, there are other aspects 
of the parameter space that must be taken into considera- 
tion, like the position of the observer, the mode degree and 
the oscillation frequency. Any definitive conclusion about a 
given observed phenomenon requires that the whole param- 
eter space is explored. Results of such an analysis are to be 
presented in paper II. 
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APPENDIX A: ORDER OF MAGNITUDE ANALY- 



SIS 



Slow Acoustic Component 



To derive equations (|24[1 and ((25}, we have considered 
the order of magnitude of each term in equations (|19l) and 
(|20p . after separating the time and the slow and fast com- 
ponents of the displacement. For that we will take as a ref- 
erence the solution that we would obtain in the absence of 
the magnetic field. The latter is of the type, 

-^i — it? ' * i.k.z 



p 



1/2' 



+ 



P 



1/2' 



(66) 



with A and B being constant amplitudes and, 



e = (67) 

where lu 2 is the acoustic cutoff frequency defined to be uj c = 
c/ (2H) and c 2 = ■yp o / po is the sound speed in an isothermal 
atmosphere iGoughl (l993l ). In the above, A — if ui < ui c 
and B — if u) > ui c - 

We may expect that the solution for the slow acoustic 
component of the displacement will vary on a scale similar 
to the solution in the absence of the magnetic field. Based 
on this assumption and considering first the case in which 
B x w B z , we can study the order of magnitude of the terms 
in equation (|20[). For that we use O to denominate the or- 
der of magnitude, and Ti to denominate the ith term in the 
equation, with i increasing from the left to the right. We 
start by noting that the non-magnetic solution varies typi- 
cally in a scale of H, i.e., 0(^ L ) « H^ 1 ^. For frequencies 
not too far from u) c , as those seen in roAp stars, this means 

that 0(^f ) » 0((u/c)M = 0(aVpo/(7Po)&0- 

Having this in mind we find that for equation ()20p the 
order of magnitude of the terms is as follows, 



Next we determine the order of magnitude of the terms 
in equation (|72[) . Also here we start by considering the case 
of B x « B z . Taking into account equation (|70p we can show 
that 0{Ts) and 0{T±) in equation ([72j) are as follows, 



0(T 3 ) 



7Po 



0(T 4 ) « O 



1 dp 
Po dz 



IPO 



IPO 



Thus, equation (I72[) can be written as equation (|24l) . 

The order of magnitude analysis presented above as- 
sumed that B x w B z - However, using the solution for £|| s 
found in this way (cf. equation (|26|> ). and returning to equa- 
tions (|19[) and (|20p . we can show that this solution is also 
appropriate in the cases of B x tending to zero and B z tend- 
ing to zero. The former is equivalent to considering the case 
in the absence of the magnetic field, so it is necessarily in 
accordance with the analysis made previously. Let us now 
consider the case of B z tending to zero. In this case the last 
term on the right hand side of equation (|27[) dominates and, 
therefore, k 2 w "^Jg-i . Going back to the analysis of the 
order of magnitude, in equation (|20p we have, 



0(2i 
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uj 2 Po \B \ 
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0(T 2 ) « O ( 

7Po 



(73) 



(74) 



Thus, we can conclude, as in the previous case, that 0(Ti) 
is always much larger than 0(T2), leading back to equation 
(|25|l . Moreover, from the latter we can conclude that in this 
case, 



O(Ti) 



0(T 2 



7Po 



( P^iu 



IPo 



(68) 



(69) 



In the magnetically dominated region, j3 « 1 and, con- 
sequently, O(Ti) >> 0(T 2 ). Thus, equation PU} leads, after 
separating the time and considering the slow component of 
the displacement, approximately to equation (|25|l . Looking 
at the latter we can conclude, under the assumptions previ- 
ously made, that 



£l.«OG8€||.)- 
Therefore, we find that 



£u « e 



(70) 



(71) 



which was to be expected, since in a magnetically dominated 
region the slow displacement is almost acoustic and, hence, 
takes place essentially along the magnetic field lines. 

Combining equation (|25[) with equation (|19[) . after sep- 
arating the time and considering the slow component alone, 
and re-organizing the terms, we find, 



d 2 £|| s dpo df|| a 



dz 2 



+ 



+ 



uj 2 p 
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Podz dz 7po S| 



dpp B x d£±, 
Podz Bz dz 
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(75) 



Bo| 2 

Therefore, we find, as before, 

£±s << £|| s . 

Let us now compare 0{Tz) with 0(T/t) in equation ([72 
when B z tends to zero, 



(76) 



0(T 3 



0(T 4 ) 



7Po 



ypo \B \ 



Therefore, also in this case, O(Ts) is larger than 0(T4), and 
we can neglect the fourth term in equation (|72l) . which leads 
us back to equation (|24|l . So, we may conclude that the ap- 
proximate solution expressed by equation (|26p. derived when 
considering B x m B z , is in fact an approximate solution for 
fii s in the more general case, since the terms that were ne- 
glected in the first case remain small, for this solution, when 
the condition B x « B z is no longer satisfied. 

Fast Magnetic Component 

To derive equations (|3ip and (|32p we again consider 
equation (|19p . after separating the time and the fast and 
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slow components of the displacement. For that we shall keep 
in mind the condition (£, S i)~ 1 d£ si /dz » (£f i )~ 1 d£fi/dz 
that was used to distinguish the above mentioned compo- 
nents. Since 0((^ a i)~ 1 d^ s i/dz) w (3(ui^J po/jpo), for equa- 
tion (|19[) we find, 

0(Ti) « 0(T S ), 
0(T 3 ) « 0(T 5 ). 

Thus, equation (|19p leads, for the fast component of the 
displacement, approximately to equation (I31|l . In addition, 
comparing the terms in the latter equation we find that, in 
all cases, the fast component of the displacement perpendic- 
ular to the magnetic field is much larger than the fast com- 
ponent of the displacement parallel to the magnetic field, 
i.e., 

Therefore, we can say that the fast displacement is al- 
most perpendicular to the direction of the magnetic field. 
This is not surprising, since in the magnetically dominated 
region this component of the displacement is expected to be 
essentially a magnetic wave, only slightly modified by the 
gas pressure. 

Finally, neglecting Ai*, when compared to £±f, we find 
that equation (|20|) leads, for the fast component, to equa- 
tion p?)) . 



